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Taking advantage of the all-sky coverage and broad frequency range of the Planck satellite, we study the Sunyaev-Zeldovich (SZ) and pressure 
profiles of 62 nearby massive clusters detected at high significance in the 14-month nominal survey. Careful reconstruction of the SZ signal 
indicates that most clusters are individually detected at least out to Rsoq. By stacking the radial profiles, we have statistically detected the radial 
' SZ signal out to 3 x Rsoo, i.e., at a density contrast of about 50-100, though the dispersion about the mean profile dominates the statistical errors 
f^*) | across the whole radial range. Our measurement is fully consistent with previous Planck results on integrated SZ fluxes, further strengthening the 
£Nj . agreement between SZ and X-ray measurements inside Rsoo- Correcting for the effects of the Planck beam, we have calculated the corresponding 
t-H ■ pressure profiles. This new constraint from SZ measurements is consistent with the X-ray constraints from XMM-Newton in the region in which 
r* 1 the profiles overlap (i.e., [0.1-1] R500), and is in fairly good agreement with theoretical predictions within the expected dispersion. At larger radii 
# I the average pressure profile is shallower than the predictions. Combining the SZ and X-ray observed profiles into a joint fit to a generalised 
k>( . pressure profile gives best-fit parameters [Po, C500, 7, c,/3] = [6.41, 1.81,0.31, 1.33,4.13]. Using a reasonable hypothesis for the gas temperature 
■ in the cluster outskirts we reconstruct from our stacked pressure profile the gas mass fraction profile out to 3i?5oo- Within the temperature driven 
(-3 uncertainties, our Planck constraints are compatible with the cosmic baryon fraction and expected gas fraction in halos. 

Key words. Cosmology: observations - Galaxies: clusters: general - Galaxies: clusters: intracluster medium - Submillimeter: general - X-rays: 
general 



1. Introduction 

In a pure hierarchical gravitational collapse scenario in the con- 
cordance ACDM cosmology, concentrations of matter ("halos") 
are fully characterised by their redshift and mean matter den- 
sity, which in turn are related to the power spectrum of initial 



* Corresponding author: E. Pointecouteau, 
etienne . pointecouteau@irap . omp . eu 



density fluctuations (Peebles 1980). The scale-free dark matter 
(DM) collapse drives the evolution of halo concentration across 
cosmic time (see e.g., recent work by Bhattacharya et al. 201 1; 
Gao et al. 2012), and the ensuing similarity yields both a univer- 
sal dark matter distribution and simple global scaling relations 
that should describe the entire halo population (Bertschinger 
1985; Kaiser et al. 1995; Navarro et al. 1995; Evrard et al. 1996; 
Navarro et al. 1997; Voit 2005; Arnaud et al. 2005). However, 
the observable properties of clusters are determined by the visi- 
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ble baryonic component, which is subject to more complex phys- 
ical processes related to galaxy formation and feedback. As the 
main baryonic reservoir in massive halos, the hot gas in the intra- 
cluster medium (ICM) is the natural target for studying the phys- 
ical processes at play and their link to the underlying cluster 
DM content. Modelling and understanding the baryon physics 
and disentangling the effect of various feedback processes is a 
very challenging task (see e.g., Borgani & Kravtsov 2009, for a 
review), underlining the need for appropriate observational con- 
straints. 

The ICM attains X-ray emitting temperatures due to grav- 
itational heating within the halo potential well. X-ray emis- 
sion is proportional to the square of the gas density, thus it 
probes the denser regions of the hot gas. The thermal Sunyaev- 
Zeldovich (SZ) effect (Sunyaev & Zeldovich 1972), due to in- 
verse Compton scattering of Cosmic Microwave Background 
(CMB) photons by the ICM, is proportional to the thermal gas 
pressure integrated along the line of sight. It is sensitive to den- 
sity and/or temperature variations, such as shocks and compres- 
sion. These two independent observational probes are thus com- 
plementary and allow us to further constrain the physics of the 
ICM. 

Recent X-ray observations based primarily on representa- 
tive samples have returned a consistent picture of the scaling 
and structural properties of halos, from high mass clusters down 
to the low mass group regime (see, e.g., Bohringer et al. 2007; 
Vikhlinin et al. 2006; Croston et al. 2008; Pratt et al. 2009; 
Arnaud et al. 2010; Sun et al. 201 1; Sun 2012). In parallel, SZ 
observations with instruments such as the Sunyaev-Zeldovich 
Array (SZA, Muchovej et al. 2007), the South Pole Telescope 
(SPT, Carlstrom et al. 201 1), the Atacama Cosmology Telescope 
(ACT, Swetz et al. 201 1), and Planck (Tauber et al. 2010) 1 have 
recently started to deliver on the promise of SZ observations for 
cluster studies. 

Building on earlier works on smaller cluster sam- 
ples (Benson et al. 2004; Andersson et al. 2011), recent re- 
sults from Planck have underlined the consistency be- 
tween the X-ray and SZ view of the ICM within Rsoq 2 . 
These constraints were achieved using three different ap- 
proaches. The first, detailed in Planck Collaboration X (2011), 
involved bin-averaging of an X-ray selected sample from 
the Meta-Catalogue of X-ray detected Clusters of galax- 
ies (MCXC, Piffaretti et al. 2011). The second, described in 
Planck Collaboration XI (2011), concerned the comparison of 
SZ measurements of 62 local clusters from the Early Release 
Compact Source Catalogue (ERCSC, Planck Collaboration VII 
201 1; Planck Collaboration VIII 201 1) with good quality XMM- 
Newton archive data. Finally Planck Collaboration IX (201 1) ex- 
amined the scaling properties of 21 newly-detected Planck clus- 
ters confirmed by XMM-Newton. These studies have provided 
well-constrained scaling relations (e.g., F500-M500, Jsoo-^x.soo) 
in the local Universe, including the first measurement of their 
intrinsic scatter, to be used as a reference for future evolu- 
tion and cosmology studies. Other dedicated investigations of 
smaller samples (e.g., Andersson et al. 201 1; Sifon et al. 2012) 

1 Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 

2 The quantity Rsao corresponds to a total density contrast 6 = 500, 
as compared to p c (z), the critical density of the Universe at the cluster 
redshift. It is linked to the mass scale by: M500 = (47r/3) 500 p c (z) R 3 5m - 



have reached conclusions similar to Planck early results regard- 
ing the agreement between the X-ray and SZ view of the ICM. 

Beyond global properties, similarity of shape in cluster ra- 
dial quantities has long been shown in X-ray observations (e.g., 
Arnaud et al. 2001; Pointecouteau et al. 2005; Vikhlinin et al. 
2006; Croston etal. 2008; Arnaud et al. 2010), in the optical 
(e.g., Rines & Diaferio 2006; Wojtak & Lokas 2010), and with 
weak-lensing (e.g., Postman et al. 201 1). For SZ studies, the ra- 
dial distribution of the thermal ICM pressure in massive halos is 
of particular interest. Recent X-ray observations of REXCESS, 
a representative sample of nearby objects, have shown that the 
ICM pressure distribution, when scaled appropriately, follows an 
approximately "universal" shape (Arnaud et al. 2010, A10 here- 
after) up to a radius of R500. The small differences in shape for 
the pressure profiles between relaxed and unrelaxed clusters, es- 
pecially in the central regions, does not seem to have an impact 
on the integrated thermal content within R500, or on other global 
properties, as noted by Planck Collaboration XI (201 1). 

A10 also compared their X-ray pressure profiles to those pre- 
dicted from numerical simulations. Although limited to 
good agreement was seen between radial pressure profile ob- 
servations and predictions from various numerical simulations 
within this radius. However, there are presently few observa- 
tional constraints beyond and consequently in this re- 
gion the shape of the "universal pressure profile" was extrap- 
olated according to the predictions from numerical simulations 
(Borgani et al. 2004; Nagai et al. 2007b; Piffaretti & Valdarnini 
2008). 

The thermodynamical state of the gas beyond /?5oo bears the 
signature of the complex physics taking place in the outer parts 
of the clusters. Characterisation of the gas in cluster outskirts 
is also necessary to unveil the level of thermal pressure in the 
cluster periphery, thus constraining their dynamical state. These 
constraints are crucial for our understanding of the formation 
and evolution of massive halos. In the cluster surroundings, ac- 
cretion along filaments is three-dimensional and non-spherical 
(Tozzi & Norman 2001) and the initial conditions of the ac- 
cretion shock are driven by the thermodynamical state of the 
(pre-shocked) in-falling material (Voit et al. 2002, 2003). The 
physical origin of any possible pre-heating is still unclear, but 
it may smooth the continuous accretion of substructures. Thus 
the entropy production at the accretion shock may be boosted 
(Voit & Ponman 2003; Borgani et al. 2005). Feedback may also 
improve the degree of thermalisation of the ICM after the shock 
(with a residual kinetic energy of about 10% beyond R500, Kay 
2004). 

Observations with X-ray telescopes have only recently 
started to provide insight into the physical properties of the 
gas in the cluster outskirts beyond ^500 (George et al. 2009; 
Reiprich et al. 2008; Urban et al. 201 1). The steep decline of the 
X-ray emission with cluster-centric radius makes such observa- 
tions extremely challenging with current instruments. A deep X- 
ray study of the Perseus cluster with Suzaku by Simionescu et al. 
(201 1 ) indicated that the gas mass fraction in the cluster outskirts 
was overestimated, possibly due to the effect of gas clumping 
and the non-virialised state of the ICM in these regions, as sug- 
gested in some numerical simulations (Nagai et al. 2007b). 

The different sensitivity of the SZ effect to the radial 
ICM distribution means that SZ observations have the po- 
tential to contribute greatly to the discussion on cluster out- 
skirts. The radial pressure distribution of the first SZ clus- 
ter samples have recently been presented based on observa- 
tions by SPT (Plagge et al. 2010), ACT (Sehgal et al. 2010) 
and SZA/CARMA (Bonamente et al. 2012). These studies con- 
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firmed that the ICM properties as seen by SZ and X-ray obser- 
vations are consistent at least out to R^m- Beyond this radius the 
SZ effect offers the interesting possibility to further constrain the 
thermal pressure support. Planck is highly competitive in this re- 
gard. It is the only SZ experiment with a full sky coverage, able 
to map even nearby clusters to their outermost radii and offering 
the possibility of an in-depth statistical study through the combi- 
nation of many observations. In this paper we present constraints 
on the thermal pressure support derived using SZ observations 
from the Planck survey. Following our previous methodology 
for scaling relations, we investigate these issues from a statistical 
point of view, working with data from 62 local clusters selected 
from the Planck ESZ sample (Planck Collaboration VIII 201 1), 
for which there are good quality XMM-Newton archival data. 

The paper is organised as follows. In the next section we describe 
the Planck mission and data, together with the sample of galaxy 
clusters used in this study. In Sect. 3 we recall the basic formal- 
ism on the SZ effect and the parametrisaton of the cluster pres- 
sure profile we use throughout the paper. Sect. 4 is devoted to a 
detailed description of the processing involved in the reconstruc- 
tion of the SZ and pressure profiles from Planck data. In Sects. 5 
and 6 we present the stacked profile of our sample and its best 
analytical representation. We discuss the comparison with other 
observational and theoretical constraints in Sect. 7, before pre- 
senting our conclusions in Sect. 8. 

Throughout the paper we adopt a ACDM cosmology with 
H = 70 km s _1 Mpc~', Q M = 0.3 and Q A = 0.7 . The quan- 
tity E(z) is the ratio of the Hubble parameter at redshift z to its 
present value, Hq, i.e., E(z) 2 = £2m(1 + z) 3 + Q\. 

2. Data 

2.1. Planck data 

Planck (Tauber et al. 2010; Planck Collaboration I 201 1) is the 
third generation space mission to measure the anisotropy of the 
CMB. Planck observes the sky in nine frequency bands covering 
30-857 GHz, with high sensitivity and angular resolution from 
31' to 5'. The Low Frequency Instrument (LFI; Mandolesi et al. 
2010; Bersanelli et al. 2010; Mennella et al. 2011) covers the 
30, 44, and 70 GHz bands with amplifiers cooled to 20 K. 
The High Frequency Instrument (HFI; Lamarre et al. 2010; 
Planck HFI Core Team 2011a) covers the 100, 143, 217, 
353, 545, and 857 GHz bands with bolometers cooled to 
0. 1 K. Early astrophysics results, based on data taken between 
13 August 2009 and 7 June 2010 (Planck HFI Core Team 
2011b; Zacchei et al. 2011), are given in Planck Collaboration 
VIII-XXVI 2011. Intermediate astrophysics results are now 
based on data taken between 13 August 2009 and 27 November 
2010. 

We used the full sky maps in the nine Planck frequency 
bands provided in HEALPix (Gorski et al. 2005) N s a e 
20- 1H full resolution. An error map associated to each fre- 
quency band is obtained from the difference of the first half 
and second half of the Planck rings for a given position of 
the satellite. The resulting jack-knife maps are basically free 
from astrophysical emission, whilst being a good represen- 
tation of the statistical instrumental noise and systematic er- 
ror. We adopted a circular Gaussian as the beam pattern for 
each frequency, as described in Planck HFI Core Team (201 lb) 
and Zacchei et al. (2011). Uncertainties in flux measurements 
due to beam corrections, map calibrations and uncertainties in 
bandpasses are expected to be small, as discussed extensively 



in Planck Collaboration VIII (2011), Planck Collaboration IX 
(201 1) and Planck Collaboration X (201 1). 

2.2. Cluster sample 

The first all-sky coverage by the Planck satellite 
(Planck Collaboration I 2011) has allowed the detection of 
dozens of clusters via their SZ signature on the CMB. The Early 
Release SZ sample (ESZ, Planck Collaboration VIII 2011) 
comprises 189 clusters of galaxies characterised over the sub- 
millimetre to centimetre wavelength range. The sample consists 
of SZ clusters and candidates detected with signal-to-noise 
ratios (S/N) spanning from 6 to 29 in the first all-sky survey. 
The sample was thoroughly validated by Planck internal quality 
assessment, external X-ray and optical data cross-correlations, 
and a multi-frequency follow-up programme for confirmation. 
The ESZ-cluster sample spans over a decade in mass, from 0.9 
to 1 5 x 1 14 M , which is essentially the full galaxy cluster mass 
range. 

In the following, we focus on a sub-sample of 62 clusters 
from the ESZ catalogue. This sub-sample is defined and ex- 
tensively described in Planck Collaboration XI (201 1), where it 
was used to calibrate scaling relations between the SZ and X- 
ray cluster properties. All 62 clusters were already known in X- 
rays (found in the MCXC, Piffaretti et al. 201 1), and have good 
XMM-Newton archive data, allowing for a high quality X-ray 
data analysis. Masses and radii for the sample were estimated 
using the M 50 q-Yx relation of A10 (see also Pratt et al. 2010), 
assuming standard evolution (see Eq. 9). 

The radius R500 was calculated iteratively as described in 
Kravtsov et al. (2006). While this sample is neither representa- 
tive nor complete, it is the largest, highest-quality SZ-X-ray data 
set currently available. The majority of objects lie at a redshift 
lower than 0.3 (and all have z < 0.5) and cover approximately 
a decade in mass. In angular size, i.e., #500, they range between 
3.7 and 22.8 arcmin, with a median value of 7.6 arcmin. These 
extended, SZ-bright clusters are thus ideal targets to investigate 
the spatial distribution of ICM thermal pressure support in clus- 
ters of galaxies by means of the spatial distribution of their SZ 
signal. 



3. Basic formalism 

3. 1 . The Sunyaev-Zeldovich effect 

The inverse Compton effect for thermalised electrons under a 
blackbody radiation field is treated in the Kompaneets (1957) 
equation. It was studied and characterised for the case of thermal 
electrons in clusters of galaxies by Sunyaev & Zeldovich (1970) 
and Sunyaev & Zeldovich (1972). The ensuing spectral distor- 
tion of the CMB spectrum in the direction of clusters is named 
the Sunyaev-Zeldovich (SZ) effect. 

In the following, we leave aside the kinetic SZ effect, a 
doppler effect, resulting from the cluster peculiar motion within 
the comoving reference frame of the Hubble flow. Also, the 
relativistic corrections on the SZ spectrum, of the order of 
kT e /m e c 2 , are not relevant for our study and are thus neglected 
(see, e.g., Pointecouteau et al. 1998; Challinor & Lasenby 1998; 
Sazonov & Sunyaev 1998). 

The intensity of the SZ effect is characterised by the dimen- 
sionless Comptonisation parameter y, the product of the average 
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fractional energy transferred per collision and the average num- 
ber of collisions: 



y 



^_ r 

m e c z J 



P(l)dl, 



(1) 



where <x x is the Thomson cross-section, m e the mass of the elec- 
tron and c the speed of light. P is the pressure produced by the 
plasma of thermal electrons along the line of sight. 

Assuming the clusters are spherically symmetric, we can ex- 
press the profile in the Comptonisation parameter as a geometri- 
cal projection of the spherical pressure profile along the axis of 
a cylinder: 



y(r) = 



cr T f Rb 2P(r')r'dr' 



m e c £ J r 



(2) 



As y is dimensionless we have y(8) = y(r). Actual SZ mea- 
surements derive from the convolution of the y-profile on the sky 
with the instrument spatial response, fpsF- 



(3) 



If the SZ brightness, which is proportional to y, is indeed 
independent of redshift, the SZ flux is proportional to the in- 
tegrated Comptonisation parameter and thus depends on the 
source distance via 



Y(0) = D A (zfY(R) = D A (zY 



CTj 

?n e c : 



f 

Jo 



2ny(r)rdr, 



(4) 



where D A is the angular diameter distance. 

In the following, the observed value of the integrated 
Comptonisation parameter is given: in units of Mpc 2 when ex- 
pressed in the source intrinsic reference frame, i.e., Y(R); and in 
arcmin 2 when expressed in the Planck satellite reference frame, 
i.e. Y(8). Analogously we can define the total integrated pressure 
within a sphere of radius r, and express it in Y units. We define 
this "pseudo" or "spherical" integrated Comptonisation parame- 
ter as 



Y sph (r) 



m e c x 



f 

Jo 



AnP(r)r 2 dr. 



(5) 



3.2. The scaled pressure profile 

As mentioned in Sect. 1, in a hierarchical scenario of structure 
formation the halo population is self-similar in scale and struc- 
ture. Profiles of physical quantities are universal once scaled ac- 
cording to their radius and reference quantities defined at a given 
density contrast, 6. Within this self-similar framework, we adopt 
throughout the paper the value of 5 — 500. The scaled pressure 
profile thus reads 



fracF(x) = 



with 



^500 



(6) 



We adopt for the pressure profile the analytical formulation 
given by Nagai et al. (2007a) for the generalised Navarro-Frenk 
White (Navarro et al. 1997; Nagai et al. 2007a, GNFW,) profile 



(*) = 



Po 



(csoo^fl + Ccsoox^F-^' 



(7) 



where x = r/R^oo, and the model is defined by the following 
parameters: Po, normalisation; C500, concentration parameter de- 
fined at the characteristic radius R500; and the slopes in the cen- 
tral (x <K I/C500), intermediate (x~ I/C500) and outer regions 
(x » I/C500), given by y, a and /3, respectively. 



A10 have fitted this analytical profile to a combination of 
observed pressure profiles derived for the REXCESS sample 
(Bohringer et al. 2007) from high quality XMM-Newton data, 
together with three sets of predicted profiles from numeri- 
cal simulations of structure formation implementing DM and 
baryon physics (i.e., radiative cooling and recipes for feedback). 
Measurements of the gas mass (integrating the density profile) 
and the X-ray spectroscopic temperature lead to the quantity 
^x,500 = A^gas,500?x,500- This parameter links to the actual in- 
tegrated parameter as seen from the SZ, i.e., F500, as 



CTj 



1 



'500 



IXSZ" 



e C z jU e '«p 



fX,500> 



(8) 



where fi e = 1.148, is the mean molecular weight of electrons 
for a 0.3 solar abundance plasma and m p is the proton mass. 
Here A X sz = 0.924 ± 0.004 from equation 19 of A10 derived 
from the REXCESS sample, and A XS z = 0.95 ± 0.04 from 
the best fit SZ scaling relation between Y500 and ix.soo from 
Planck Collaboration XI (201 1) using the same sample as in this 
work. For consistency purposes we use this latest value in the 
following. The predicted and measured values for Igoo are com- 
patible within the ±lcr limits. 

ix.soo is a good proxy for the cluster total mass 
(Kravtsov et al. 2006; Nagai et al. 2007a; Arnaud et al. 2007). 
The non-standard scaling relation fitted against the REXCESS 
data is provided by their equation 2: 



E{zf l5 M 5m = 

jq!4.567±0.010 



2x 10 14 M o keV 



0.561±0.018 



Mr. 



(9) 



In turn the characteristic pressure P500 scales with the cluster 
total mass, reflecting the mass variation expected in the standard 
self-similar model, purely based on gravitation (Eq. 5 in A10): 



500 



= 1.65 X 10" 3 h(zf> ' 3 
M500 



3 x 10 14 /i-! M, 



70 



!/3 



htjgke V cm 3 . 



(10) 



The M-Yx relation given in Eq. 9 deviates from the standard 
self-similar case (which has a slope of 3/5). From the definition 
of P500, any deviation from the standard self-similar scaling will 
appear as a variation of the scaled pressure profiles. As shown 
by A10 in their equation 9, this variation can be expressed as a 
function of the total mass. At S = 500, this mass dependence is 
almost constant with radius, and it can be approximated by 



P(r) 
P500 



= F(x) 



M500 



3 x 10 14 /i;'M c 



0.12 



(11) 



3.3. Conversions and normalisations 



In the following, we measure SZ profiles from the Planck data 
(Sect. 4.2). All through this work, we compare these observed 
profiles with the predictions that arise from the X-ray con- 
straints. For all clusters we assume the universal pressure profile 
shape as derived from A 10. We parametrise it according to the 
quantities derived from the X-ray analysis R^m and ix.soo (see 
previous section). We can then derive predicted pressure pro- 
files (see Eq. 7), and Comptonisation parameter profiles (Eqs. 2 
and 4). The later convolved by the instrument beam (Eq. 3) is 
directly comparable to our observed Planck profiles. 
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We apply a statistical approach to the observed SZ profile 
and to the pressure profile, averaging individual profiles in our 
sample once scaled. Profiles are scaled in radius according to 
Rsoo- The observed SZ profiles (or predicted SZ profiles from X- 
ray constraints) are dimensionless, and we therefore normalise 
them by the quantity 

<S> = Y 500 /R 2 m . (12) 

In the following we used the value of Ysqq as defined in Eq. 5. 

To translate our measured SZ profiles, we deconvolve and 
deproject them for each cluster following Eqs. 2 and 3. In prac- 
tice this operation is performed as described in Sect. 6.1. The 
resulting profiles are converted to pressure such that 

m P c 2 1 d6 
P{r) - y{6)-. (13) 

<x T D A (z) dr 



unavoidably leads to an oversampling of the pixels or, more pre- 
cisely, to a redundancy of the original all-sky map pixels in the 
reprojected patch. We produced an associated patch that tracks 
this redundancy. We also extracted equivalent patches from the 
associated all-sky SZ error map, and from the all-sky variance 
map. 

We then computed a profile from each SZ map patch. These 
were calculated on a regular radial grid with bins of width 
Ar//?5oo = 0.25, allowing us to sample each cluster profile with 
four points within Rsoo- The y value of a bin was defined as the 
mean of the values of the pixels falling in each annulus. We sub- 
tracted a background offset from the maps prior to the profile 
computation. The offset value was estimated from the surround- 
ing region of each cluster where r > 7 ^500 ■ The uncertainty 
associated with this baseline offset subtraction was propagated 
into the uncertainty of each bin of the radial SZ profile. 



4. Reconstruction of the SZ profile 

4.1. SZ signal extraction methods 

The thermal SZ maps were recovered from a combination of 
Planck channels, making use of either real or Fourier space 
methods. We applied three methods based on internal lin- 
ear combination (ILC) algorithms to the Planck data: (1) the 
Modified Internal Linear Combination Algorithm (MILCA, 
Hurieretal. 2010);, (2) Needlet Internal Linear Combination 
(NILC, Delabrouille et al. 2009); and (3) the Generalized 
Morphological Component Analysis (GMCA, Bobin et al. 
2008). A description of each method is given in App. A. 

To optimise the reconstructed SZ map we have to take into 
account that: (i) the brightness of the SZ effect (as an incre- 
ment or decrement) is maximum in the sub-millimetre to mil- 
limetre range; (ii) the typical angular size of the cluster is a 
few arcmin; and (iii) the final resolution of the reconstructed 
SZ map is determined by the lowest resolution of the combined 
frequency maps. We therefore restricted ourselves to the use of 
the six Planck-HFl channels in the SZ map reconstruction pro- 
cess. The final resolution of the reconstructed SZ maps is that of 
the lowest-resolution 100 GHz channel, or 10 arcmin. SZ maps 
are in units of the (dimensionless) Comptonisation parameter. 
Integrated Comptonisation parameters are expressed in arcmin 2 
for the observed values and in Mpc 2 for the intrinsic value in 
the source reference frame. The associated SZ noise maps are 
built in a fully consistent way from the individual frequency er- 
ror maps (see Sect. 2). 

4.2. SZ profile computation 
4.2.1. Profile extraction 

Individual cluster profiles were computed from the reconstructed 
all-sky SZ maps. Uncertainties were obtained from the all-sky 
reconstructed SZ error maps, which were derived by applying 
the same reconstruction methods to the frequency error maps 
(see Sect. 2). For instance, in the case of a linear method (e.g., 
ILC based), the coefficients used for the linear combination of 
the frequency maps were propagated to the combination of the 
frequency error maps. 

We extracted a square patch of side 20 x 6500 around each 
cluster position from the all-sky SZ map. Patches were projected 
from HEALPix (Gorski et al. 2005) to a tangential projection. 
For each patch, the pixel size was adapted so that it was con- 
stant in scaled units of 6500 over the full cluster sample. This 



4.3. Expected SZ profile 

We used the GNFW pressure profile shape characterised by A10 
from the REXCESS sample and numerical simulations to gener- 
ate a pressure model for each cluster of our sample. We kept the 
A10 best fit values of C500 = 1-18 and the three slopes, a = 1.05, 
/? = 5.49, y = 0.31. We derived the normalisation Pq from the 
observed value of ix.soo (see Eq. 5), using a conversion of 0.95 
between Fx and Fsz as described in Sect. 3.2. Finally ^500 was 
fixed to the values derived from the XMM-Newton analysis. Both 
Yx and ^500 are reported in Planck Collaboration XI (201 1). We 
computed the projection matrix and the PSF redistribution ma- 
trix (for a Gaussian beam of 10 arcmin FWHM) as expressed 
in Eqs. 2 and 3, respectively, to multiply each pressure profile 
model. The derived y profile model is directly comparable to 
the Planck measurements. The average profile model across the 
sample was derived similarly to the observed stacked profile (see 
above), and is used to compare to the observed stacked profile. 

4.3.1. Accounting for correlation between points 

A certain level of correlation is introduced between the points of 
a radial profile derived in this manner. We account for this cor- 
relation in the covariance matrix of each profile, which is com- 
puted as follows. 

For each map patch, we masked radii 9 < 7 #500 centred on 
the cluster and computed the power spectrum of the noise on the 
remaining area. The use of the surrounding regions of the cluster 
to characterise the noise properties allows us to account for the 
effect of astrophysical contamination in the cluster vicinity, as 
well as the instrument noise and systematics. These sources of 
contamination were, by construction, excluded from the all-sky 
frequency error maps (see Sect. 2). The power spectrum drawn 
from the SZ error map of each patch is therefore, systematically, 
slightly lower than the one computed over the cluster surround- 
ings. 

We then simulated m = 500 realisations of the noise patch 
based on this power spectrum, accounting for the variance map 
and sky pixel redundancy and assuming an inhomogeneous cor- 
related Gaussian noise. We then extracted a noise profile from 
each realisation, reproducing the baseline background subtrac- 
tion used for the observed profile. The covariance matrix was 
built from all the simulated noise profiles, i.e., C = Pj,P„, 
where P„ is an n points x m matrix of simulated noise profiles. 
Considering two points as correlated when their correlation co- 
efficient is larger than 0.3, the typical level of point-to-point cor- 
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Fig. 1: Average scaled SZ profiles for each of the nine Planck frequency bands (decreasing frequency from left to right and from top to bottom). 
The points within each individual profile are correlated at about the 20% level (see Sect. 4.3.1). At each frequency, the stacked radial profile 
is obtained from the average of the 62 individual cluster profiles scaled by R500 and 0,/(0,), in the x- and v-axis directions, respectively (see 
Sect. 5.1). The light-red shaded area in each panel indicates the statistical uncertainty about the average observed profile. The solid black line 
corresponds to the prediction of the universal pressure profile (A10) for our cluster sample in each band, and the two dashed lines indicate its 
associated dispersion. 



relation in our profiles is about 20% (16%, 21% and 28% for 
MILCA, NILC and GMCA, respectively). 

4.3.2. Stacking procedure 

Given that all profiles are computed on the same grid in scaled 
radius, they can easily be stacked together. Each profile y,- was 
rescaled by R^xi.i and <t>, (as given in Eq. 12), respectively in 
the X and Y axis directions. The stacked profile and associated 
covariance matrix are given by 

y^Vj and C=if% (14) 
n L-i O; n 2 (J) 2 

i it 



For the computation of <t>, we used Y sp h$M in Mpc 2 as given 
in Eq. 5. For simulated data (App. B), it was derived from the 
Planck M500-F500 relation (Planck Collaboration IX 2011); for 
observed data (Sect. 5.3), it was measured directly from the 
Planck data. 

When for a given bin, values for all clusters are dominated 
by the signal, we assumed a log-normal distribution of their scat- 
ter and stacked them in logarithmic space. The logarithm of the 
stacked value and associated error were then translated back into 
linear space. In other terms, for bin j, if all the clusters satisfied 
the condition (yj - cr J y ; ) > 0, we stacked all the measurements in 
the logarithmic space; otherwise, we did it in the linear space. 
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We computed the statistical and weighted (i.e., 1/cr 2 ) aver- 
age profiles, as well as the median profile of the sample, and 
checked their consistency. Over a radial range of 3 xR^qo, the av- 
erage of the relative error of the weighted mean and median pro- 
files with respect to the statistical average profile is smaller than 
10% and 5%, respectively. Based on this agreement, we used the 
statistical average to compute the stacked profile throughout this 
study. 

Of the 62 clusters in our ESZ-XMM sample, only two are 
spatially close and thus are potentially physically connected: 
A3528 and A3532. These two clusters are members of the 
Shapley supercluster. Each of those clusters was masked over an 
area of radius 3 R500 when processing the other. Thus, we con- 
sider that their profiles and covariance matrices are independent. 

At the end of our data processing, the stacked covariance ma- 
trix encompasses the statistical errors due to instrumental noise, 
astrophysical fluctuations at the cluster location and systematic 
effects (e.g., instrumental and arising from data processing). For 
the purposes of our study we have also computed the dispersion 
across our sample for each position within our profiles. Both un- 
certainties are propagated throughout the analysis, and are cross- 
compared in the following sections. 

5. Planck galaxy cluster SZ profile 

5.1. Frequency stacked profiles 

We first looked for the raw cluster signature in each of the nine 
Planck frequencies. We extracted a brightness profile at each fre- 
quency from the raw Planck maps for each of the 62 clusters in 
the ESZ-XMM sample. These were rescaled in units of R^o in 
the jt-axis direction and in units of <!>,■/ (<!>,■) in the y-axis direction 
(see Sect. 3.3), and then stacked. The resulting average bright- 
ness profiles are shown in Fig. 1, where the shaded area in each 
panel depicts the statistical error about the mean flux value in 
each bin. 

The cluster signal appears clearly in most of the nine fre- 
quencies, following the SZ thermal spectral signature. The sig- 
nal is positive at 545 and 353 GHz, compatible with zero at 
217 GHz, and then negative down to 30 GHz. The 62 clusters 
of our sample are strong SZ sources, as expected for objects de- 
tected at S/N > 6 (Planck Collaboration VIII 201 1). The stack- 
ing procedure mostly averages out the effect of foreground and 
background contaminants; however, they contribute to the final 
dispersion in the stacked profiles. Owing to their Gaussian na- 
ture, CMB fluctuations are more easily washed out in the stack- 
ing procedure than the Galactic and extragalactic dust emission 
in the high frequency range, and Galactic and point source radio 
emission (i.e., free-free and synchrotron) in the lower frequency 
range. 

We also plotted for each frequency the expected profile com- 
puted as described in Sect. 4.3. This comparison shows a first 
order agreement between the Planck measurements and the ex- 
pected SZ profiles derived from X-ray constraints assuming a 
GNFW pressure profile shape. 

5.2. Stacked SZ image 

Figure 2 presents the stacked average image of the 62 clusters of 
our ESZ-XMM sample, plotted on a logarithmic scale. Before 
averaging, each individual SZ map was normalised identically to 
the profiles (seeEq. 12) and randomly rotated by 0°, 90°, 180° or 
270°. The stacked image was then renormalised by (O;) in order 
to be expressed in y units. The image shows that the SZ signal is 




Fig. 2: Stacked scaled image (size is 10 ^500 on a side) using a logarith- 
mic stretch in y. Individual maps are rescaled by C>, before averaging, 
and then multiplied by <0,). The black and white circles mark the loci 
of 1 and 3 x R500, respectively. 



detected out to ~ 3 x R500, marked by the white circle in the fig- 
ure. The extent of existing observational constraints on the gas 
distribution in clusters in SZ and X-rays, i.e., R500, is marked by 
the black circle. We also generated a jackknife-type map by aver- 
aging the scaled individual SZ maps, m,-, as £ (-l)'m,-. The rms 
values of both the stacked SZ and on-off maps beyond 5 x R500 
are compatible: 4.6 x 10~ 7 and 4.8 x 10~ 7 , and, in turn, compati- 
ble with the rms value below this radius on the jackknife stacked 
map: 3.5 x 10~ 7 . This check demonstrates that our average map 
is not strongly affected by residuals or biases. The radial pro- 
file computed from the stacked map is fully compatible with the 
stacked profile discussed below. 



5.3. Observed stacked SZ profile 

We derived the observed stacked Planck SZ profiles for the 
ESZ-XMM sample as described in Sect. 4.2. T500 is obtained 
from the algorithms used for blind detection of SZ clusters in 
the Planck survey (Planck Collaboration VIII 201 1), namely the 
Powell Snakes (PwS, Carvalho et al. 2009, 2011) and multi- 
frequency matched filter (MMF3, Melin et al. 2006) algorithms. 
In both cases the algorithms implement a universal pressure pro- 
file shape (A10) with the position of the cluster fixed to the 
XMM-Newton coordinates and the size of the universal pressure 
profile fixed to /?5oo- The fluxes from both methods (i.e., PwS 
and MMF3) are consistent over the whole sample. The median 
value for the ratio PwS/MMF3 is 0.96 ± 0.05 (see App. B.2). To 
further validate the above, and as a consistency check with pre- 
vious Planck results, we fitted each individual SZ profile with 
a projected, PSF-convolved universal pressure profile. We fixed 
R500 to the best fitting X-ray value from Planck Collaboration XI 
(2011) and only fitted the normalisation, Fsoo- Given the agree- 
ment, we used the MMF3 values to compute <t> for each cluster 
(Eq. 12). 
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Fig. 3: Left: individual profiles (grey lines) and scaled stacked profile (red points) for our sample of 62 clusters. The light-red shaded area marks 
the dispersion about the average stacked profile, with its upper and lower limits highlighted by the two dashed red lines. The points within each 
individual profile are correlated at about the 20% level (see Sect. 4.3.1). The solid black line (labeled in the legend as "Expected from Paio") is 
the average stacked profile obtained from the expected SZ individual profiles drawn from the universal pressure profile of A10, parameterised 
according R^oo and J500 derived from the XMM-Newton data analysis presented in Planck Collaboration XI (2011). The two dotted black lines 
indicate the associated dispersion about this model profile. The bottom panel shows the value of x at each point with respect to the expectations 
from the universal profile taking into account the measured errors. Right: matrix of the correlation coefficient for the observed radial stacked 
profiles. Both panels are illustrated here for the MILCA results. 



Appendix B presents a detailed investigation of the conver- 
gence between stacked profiles derived from each of the three 
different methods of SZ signal reconstruction. We find that all 
three methods agree remarkably well over the entire radial range, 
both for simulated and observed SZ profiles. The three SZ re- 
construction methods lead to profiles fully compatible with each 
other. Across the range of radii over which the profiles have been 
computed (i.e., [0-10]x/?5oo), accounting for the correlated er- 
rors of each profile, the reduced x 2 of NILC and GMCA with 
respect to MILCA are 0.48 and 0.62, respectively 3 . 

The left panel of Fig. 3 shows the stacked MILCA profile 
for the 62 clusters in the ESZ-XMM sample compared to the 
individual profiles (see Sect. 4.2). The SZ signal is statistically 
detected over more than two decades in intensity, and out to a 
remarkably large radius of ~ 3 R500, reaching far into the clus- 
ter outskirts. Assuming that the virial theorem can be applied 
(which at such large radius might be breaking down), the outer 
radius of our statistical detection corresponds to a density con- 
trast of 6 ~ 50. More conservatively, we can assume that we are 
statistically probing the average SZ and pressure distribution 
down to regions of density contrast of 5 ~ 50-100. The disper- 
sion about the mean profile dominates the statistical uncertain- 
ties. It is minimal (by construction) at ~ 20% around /?5oo, but in- 
creases towards the centre and the outskirts to ~ 50% and ~ 65% 
at 0.3 and 2x/?5oo, respectively. At the most external radius of our 
statistical detection, i.e., ~ 3 xR^m, the dispersion of the individ- 



3 As our three reconstruction methods are remarkably consistent, for 
the clarity of display and discussion, hereafter we illustrate our pre- 
sentation with one of the three methods only. In each case the specific 
method used will be indicated. 



ual profiles about the mean is more than 100%, as at these large 
radii, the individual SZ profiles are fully dominated by noise. 

The right panel of Fig. 3 shows the correlation coefficient 
matrix for the stacked profiles presented on the left panel, and il- 
lustrates the degree of correlation between the points in the pro- 
files. 



5.4. Comparison with expectations from universal profile 

The stacked model used for comparison was computed as de- 
scribed in Sect. 4.3. When considering the statistical errors only, 
the measured stacked profile is significantly above the model. 
Taking into account the error on the model and the correlated er- 
rors between points of the reconstructed SZ profile, we obtained 
a reduced^ 2 value of 3.53 within 3/?5oo- If we omit the error 
on the model, this value becomes 7.48 As our tests on simula- 
tions show a very good agreement between the input model and 
the output reconstructed profiles (see App. B), the difference ob- 
served here between the measured SZ profiles and the predicted 
model from X-ray constraints is not an artefact of the method, 
but an intrinsic difference. This difference is significant at a 2- 
3cr level from around Rsqq out to 3 ^500- 



6. The galaxy cluster pressure profile from Planck 

6. 1. PSF deconvolution and deprojection 

The pressure profile is derived for each cluster in our sample by 
applying a deconvolution and deprojection to the observed SZ 
profiles. As already mentioned, we assume spherical symmetry. 
The convolution by the instrumental beam and the geometrical 
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Fig. 4: Left: Planck pressure profile obtained from the average of the individual pressure profiles across our sample (red points) shown together 
with the stacked pressure profile derived from the XMM-Newton data for the same sample (purple points). The dispersions about the SZ and X-ray 
profiles are depicted, respectively, by the red and purple shaded areas. Our best fit GNFW profile is shown as a solid black line and that of A10 
as a solid green line. The lower panel shows the x profile of these two best models taking into account the statistical errors and dispersion about 
the observed profile. Right: stacked profile from Planck and XMM-Newton for the sub-samples of cool-core (open circles) and non cool-core 
(open stars) clusters within the ESZ-XMM sample. The shaded areas are identical to the one shown in the left panel. Our best fit models for each 
sub-sample are shown as black solid and dashed lines for the cool-core and non cool-core clusters, respectively (see Table 1). The best A10 fit for 
cool-core and non cool-core clusters are shown as green solid and dashed lines, respectively. The lower panel gives the x profiles of our best fit 
models with respect to the measured profiles and associated errors (including dispersion). On both panels the data points are correlated at about 
the 20% level (see Sect. 4.3.1) 



projection of a spherical pressure profile into a Comptonisation 
parameter profile are expressed in Eqs. 3 and 2, respectively. 
We have applied a deconvolution plus deprojection algorithm 
adapted from the method described by Croston et al. (2006). 
This method allows us to perform a straightforward real space 
deconvolution and deprojection using a regularisation procedure 
originally applied to X-ray surface brightness profiles to derive 
the underlying density profile. The method is adapted to lower 
statistics (i.e., lower number of bins per radial profile). In the 
present case, the PSF response matrix was computed for a cir- 
cular Gaussian beam with FWHM = 10 arcmin, i.e., the angular 
resolution of our SZ maps. SZ cylindrical profiles deconvolved 
from the PSF were then deprojected into 3D profiles (assuming 
spherical symmetry for the clusters), and converted to pressure 
following Eq. 13. 

We used a Monte Carlo (MC) method to propagate the mea- 
sured error over the SZ profiles onto the PSF-corrected and de- 
projected profiles, accounting for the covariance between points 
expressed in the form of the matrix C, for the cluster ;. For each 
cluster i, its covariance matrix was Cholesky decomposed (i.e., 
C, = LjLj). We assumed correlated Gaussian noise, and gener- 
ated m = 10,000 realisations of the SZ profile, P,, with respect 
to this decomposition (i.e., P' { = Pi + xL{, where x is an array 
of random numbers following a normal distribution). With this 
method we make sure we properly sample the noise properties of 
each profile. Each realisation was then rebinned logarithmically 
in radius and run through our deconvolution and deprojection 
algorithm. We conservatively chose to feed to the algorithm val- 
ues of the weight per point equal to the inverse of the square of 
the flux dispersion in this bin (i.e., w = l/o~g ux ). The weights 



remained the same for all realisations. These m realisations of 
the pressure profile were used to compute the covariance ma- 
trix of the cluster pressure profile similarly to what is described 
in Sect. 4.3.1. With this MC process, we ensure a proper prop- 
agation of the correlated noise through the PSF-correction and 
deprojection of the SZ profile. 

6.2. The stacked pressure profile 

As we did for the SZ profiles, the 62 pressure profiles were 
rescaled in the x- and y-axis directions according to P500 and 
P500, and then stacked together. Following Sun et al. (201 1) we 
multiplied individual profiles by the factor f(M) = (M500/3 X 
10 14 hzl M ) 012 to further correct for the slight dependency of 
the pressure profiles with mass (Eq. 11). For our ESZ-XMM 
sample, the average of this factor is (f{M)) = 1.09. Our three 
reconstruction methods lead to compatible stacked pressure pro- 
files. However, we have accounted for their fluctuations by 
adding to the diagonal terms of the covariance matrix of the 
stacked profile the maximum point-to-point difference between 
the MILCA profile and the two other methods. In parallel to our 
Planck SZ analysis, we have also derived the pressure profiles 
of our 62 clusters from the analysis of XMM-Newton data. We 
followed the method applied to the REXCESS sample and pre- 
sented in A 10. Details of the XMM-Newton analysis and results 
are provided in Democles et al. (2012). 

The Planck and XMM-Newton stacked pressure profiles de- 
rived from our sample of 62 ESZ-XMM clusters, are displayed 
in Fig. 4. The two sets of data agree remarkably well. They 
are fully compatible within their respective dispersions, and 



9 



Planck Collaboration: Pressure profiles of galaxy clusters from the Sunyaev-Zeldovich effect 



they overlap in the radial range ~ (0.3-1) Rsqq. We recall that 
the points from the Planck data suffer a degree of correlation 
at about the 20% level, whereas all the points in the XMM- 
Newton profile can be considered as uncorrected. These two 
independently derived pressure profiles are fully complemen- 
tary. Indeed, in the central region the pressure profiles are very 
well constrained from the X-ray data, whilst the Planck measure- 
ments hardly reach down to 0.1 /?5oo due to the moderate spatial 
resolution of our SZ maps. At larger radii, the extent of the X- 
ray observations is limited to radii smaller than ~ (0.7- 1.0) ^500 
because of the XMM-Newton field of view and of the quick drop 
in X-ray surface brightness with increasing radius. The Planck 
profile extends far beyond this radius. 

With this joint constraint, we bring for the first time a com- 
prehensive observational view of the distribution of the average 
thermal pressure distribution in clusters of galaxies out to a den- 
sity contrast of 6~ 50-100. 

Following A10, we have investigated two sub-samples of 
our ESZ-XMM clusters. Keeping the cool-core (CC) ver- 
sus non cool-core (non-CC) classification as provided in 
Planck Collaboration XI (201 1), we computed the stacked pro- 
files for the 22 CC versus 40 non-CC clusters. Both stacked 
profiles are displayed in the right panel of Fig. 4. They are, as 
expected, different in the central parts, with a more peaked pro- 
file for CC systems and a shallower one for non-CC clusters. 
However, in the outer parts, our observed profiles for the two 
subsamples have very similar slopes. 

6.3. Best-fit to the Generalised NFW profile 

We have combined the Planck and XMM-Newton pressure pro- 
files into a joint fit to the generalised NFW model (Eq. 7). The 
pressure profiles derived over our sample from the Planck and 
XMM-Newton data cover a radial range from 0.02 x R500 to 
3x/?5oo, with an overlap in the radial interval of (0.1-l)x/?5oo- In 
order to find the best parametric representation of our observed 
stacked profile, we ran a Monte Carlo Markov Chain analysis to 
find the maximum likelihood solution parametrising the GNFW 
profile. The Planck and XMM-Newton data are two different and 
fully independent observational data sets. We therefore com- 
puted the likelihood as the product of the two independent like- 
lihoods: X. = -Cx x -£sz- For the Planck data, we have accounted 
for the correlation between points through the use of the covari- 
ance matrix (see Sect. 6. 1). We also accounted for the dispersion 
about the average X-ray and SZ profiles across the sample by 
quadratically adding it to the errors bars of the XMM-Newton 
points and to the diagonal elements of the covariance matrix for 
the Planck profile. 

We performed our MCMC analysis in log-log space for var- 
ious combinations of free parameters. Fixed parameters were as- 
signed the values from the best fit provided by A 10. Our best fit 
parameters are reported in Table 1 . In order to find the best possi- 
ble analytical representation of our data, we assessed the quality 
of the fit via the value of the reduced^ 2 using only the statistical 
errors, x^ m . The configuration with four and five free parameters 
are of equivalent quality, the former being slightly better (i.e., 
A^tat = ar, d 1-0, respectively). Unsurprisingly the degeneracy 
between free parameters increases with their number, translating 
into a drastic increase in the parameter covariance. We therefore 
adopted the four free parameter approach as our best parametri- 
sation of a GNFW profile of our data, i.e., [Pq, C500, a, ft] free, 
with y fixed to 0.31. The associated values are marked in bold 
face in Table 1 , and are given purposefully with high precision 
with respect to the errors on parameters. These uncertainties are 



Table 1 : Best fit parameters for a generalised NFW pressure pro- 
file, with different numbers of fixed parameters. The^- 2 tat value is 
the reduced^ 2 computed with respect to the statistical error bars 
(i.e., not including the dispersion across the sample). The bold- 
face line indicates the best parametric representation of our com- 
bined XMM-Newton and Planck pressure profile (see Sect. 6.3 
for further details). The number of degrees of freedom is given 
in the last column. 



A) c 500 


7 


a 


P 


Xstat 


d.o.f 


All clusters 


6.32 1.02 
6.82 1.13 
6.41 1.81 

5.78 1.84 


0.31 + 
0.31+ 
0.31+ 
0.35 


1.05+ 
1.05+ 
1.33 
1.39 


5.49+ 
5.17 
4.13 
4.05 


3.8 
5.8 
0.9 
1.0 


15 
14 
13 
13 


Cool-core clusters 


11.82 0.60 


0.31+ 


0.76 


6.58 


1.1 




Non cool-core clusters 


4.72 2.19 


0.31+ 


1.82 


3.62 


1.2 





+ Fixed parameters which are assigned the best fit values of the A10 
profile. 

illustrated by the corresponding marginalised posterior likeli- 
hood distributions that are shown for individual (ID) and pairs 
of parameters (2D) in Fig. 5. 

In parallel, following exactly the same procedure, we also 
fitted the stacked profiles for the sub-samples of cool-core and 
non cool-core clusters. The best fit parameters are also reported 
in Table 1 . Both the stacked profiles and best fit model are shown 
on the right panel of Fig. 4. 

We further adopted this four free parameter configuration to 
fit each of our measured individual pressure profiles, combin- 
ing XMM-Newton and Planck pressure data. As describe in the 
above section for the stacked profile, we accounted for the point- 
to-point differences between our three reconstruction methods 
(see Appendix A) in the error budget. Results are reported in 
App. C, Fig. C.l and Table C.l. 

7. Discussion 

7. 1 . The observed pressure profile 
7.1 .1 . The core of clusters 

Our present results are in very good agreement with the uni- 
versal profile derived by A10 outside the core and down to 
Rsoo- Within the core, i.e., R < 0.15/?5oo, our observed pro- 
file lies significantly below the A10 profile. In comparison to 
the REXCESS sample which is X-ray selected, our ESZ-XMM 
sample is SZ selected (except for the intersection with the XMM- 
Newton archives) and is thus, a priori, closer to being a mass- 
selected sample. It therefore contains more dynamically dis- 
turbed clusters (e.g., A2163, RXJ2228+2036, etc.). As has also 
been indicated by results in the validation follow-up of Planck 
clusters with XMM-Newton (Planck Collaboration IX 2011; 
Planck Collaboration I 2012; Planck Collaboration IV 2012), 
this seems to indicate that X-ray selection under-samples the 
population of morphologically disturbed clusters. When split- 
ting our sample into CC and non-CC clusters, the average pro- 
file for the CC clusters (and its associated best fit) is in better 
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Fig. 5: Marginalised posterior likelihood distribution for our best MCMC fit to a GNFW pressure profile with free parameters [Pq, csoo, o,P\ 
and y = 0.31. The best fit values of the parameters are marked as white crosses and solid red lines, respectively, on the 2D and ID posterior 
distributions. On the 2D distributions, the white triangles mark the value of the A10 parameters, while the red, purple and green solid lines define, 
respectively, the 68, 95 and 99% confidence levels. 



agreement with the universal profile given in A 10. However, it 
is still slightly lower in the central parts than the A10 profile for 
CC clusters (see right panel of Fig. 4). There might thus still be 
differences between the population sampled by our CC clusters 
and those from A 10. On the contrary, our non-CC cluster pro- 
file agrees well with the A10 non-CC profile, except for the very 
central parts (R < 0.04 R 50 o)- 

7.1 .2. The inner profile 

Observationally, it is hard to accurately determine the distribu- 
tion of the SZ signal and, thus, of the underlying thermal pres- 
sure. Early single target SZ studies were limited to the clusters' 
inner regions (i.e., within R500), and were not really competitive 



with the X-ray measurements in terms of resolution and sensi- 
tivity (see e.g., Pointecouteau et al. 2001; Halverson et al. 2009; 
Korngut et al. 2011). However, the combination of SZ and X- 
ray tracers has already shown its potential in terms of structural 
studies of clusters (Pointecouteau et al. 2002; Kitayama et al. 
2004; Jia et al. 2008; Basu et al. 2010). The first studies work- 
ing on SZ cluster samples were presented by the SPT and ACT 
collaborations (Plagge et al. 2010; Sehgal et al. 2010), with 15 
and 9 high significance clusters, respectively. The SPT team 
found good agreement between the X-ray predicted signal and 
their SZ measurement within /?5oo- Similarly, in a recent study, 
Bonamente et al. (2012) used a sample of 25 clusters observed 
with SZA/CARMA (Mohr et al. 2002) and found an excellent 
agreement when modelling the SZ emission over their sample, 
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either with the universal pressure profile from A10 or with the 
model from Bulbul et al. (2010). Our SZ and X-ray data are in 
excellent agreement over the radial range (0.1-1) R500, providing 
tighter observational constraints. 

Our results, together with the SPT, and ACT results, have to 
be compared with other works that claim a difference between 
SZ and X-ray measurements. It was first reported by Lieu et al. 
(2006) that the SZ effect only accounted for one fourth of the 
expected signal from X-ray constraints on the basis of a WMAP 
(one year) analysis of 31 randomly selected nearby clusters. 
A statistical analysis based on WMAP-1 data of 49 clusters 
(Komatsu et al. 201 1) found the ratio between the measured and 
expected SZ signal to be 0.5-0.7. However, these claims did 
not agree with other independent WMAP analyses (Melin et al. 
201 1). Neither do they agree with strong observational evidence 
from the Planck frequencies common to WMAP (see profiles for 
30, 44 and 70 GHz in Fig. 1), nor with the excellent match we 
find here between X-ray and SZ properties within R^m- 

7.1 .3. The outskirts 

Beyond R500, the observational situation is even more wanting. 
Only a few X-ray observations with XMM-Newton, Chandra 
and/or Suzaku constrain the density, temperature or gas frac- 
tion profiles out to ~^200 - 1-4 ^500 (George et al. 2009; 
Reiprich et al. 2008; Urban et al. 201 1; Simionescu et al. 201 1; 
Walker et al. 2012a,b). This type of X-ray detection remains 
very challenging and requires very long exposure times, as at 
larger radii the X-ray emission is extremely faint. We recall 
that beyond R500 the universal pressure profile from A10 was 
constrained by predictions from numerical simulations, not by 
observations. Probing the gas with the SZ effect is therefore a 
powerful alternative, as shown, for instance, by the SPT aver- 
age emission over 15 clusters (Plagge et al. 2010) or the WMAP 
statistical analysis over ~ 700 clusters by Atrio-Barandela et al. 
(2008), though the latter is affected by the limited resolution and 
sensitivity of the survey. Conversely, from the tentative analysis 
of WMAP-3 data for 193 clusters with X-ray temperature above 
3keV, Afshordi et al. (2007) provided constraints for the clus- 
ter pressure profile (significant out to ~ 1.3 /^oo)- In agreement 
with these earlier works, the Planck measurements of the present 
study are the first to allow a precise description of the thermal 
pressure distribution out to the cluster outskirts. 

As shown in App. C, Planck also resolves some in- 
dividual pressure profiles. The work presented by the 
Planck Collaboration V (2012) on the Coma cluster is even more 
striking, as the SZ profile detection reaches beyond 3 x /?5oo, i.e., 
as far as the statistical measurement presented in this work. The 
derived best fit model for the Planck observations of Coma is in 
full agreement with our stacked result. 

We have also investigated the change in integrated 
Comptonisation parameter due to the difference in pressure pro- 
file with respect to A10. Within a fixed aperture R500, known 
from ancillary data, the difference is marginal as the two pro- 
files are alike. The average ratio across the sample between 
the value of Y500 from our best profile and that of A10 is 
1.02 ± 0.03, consistent with the value of 1.01 derived for 
a GNFW model with the two parametrisations. This consis- 
tency demonstrates the robustness of the previously-published 
Planck SZ analysis and scaling relations where the A10 pro- 
file was adopted as a fiducial model for known Rsoa values 
(Planck Collaboration IX 2011; Planck Collaboration X 2011; 
Planck Collaboration XI 2011; Planck Collaboration I 2012; 
Planck Collaboration III 2012; Planck Collaboration IV 2012). 



The ratio for the values of Y for the two profiles is 1.13 and 1.19 
within fixed apertures of 3 R500 and 5 R500 respectively, as ex- 
pected given the slightly flatter outer slope derived in this work 
with respect to that of A 10. 

7.2. Comparison with theoretical predictions 
7.2.1. Sets of numerical data 

We have investigated three sets of simulated clusters in order 
to compare to our combined Planck and XMM-Newton pressure 
profile. All are taken from numerical simulations of structure 
formation in a ACDM cosmology: 

1. We first compared to the set of combined simulations used 
by A10 to derive their universal pressure profile together 
with the REXCESS data. They comprise 93, 88, and 14 
simulated clusters with M > 1O 14 M , from Borganietal. 
(2004), Nagaietal. (2007a) and Piffaretti & Valdarnini 
(2008), respectively. We refer to this set of simulations as 
B04+N07+P08 hereafter. 

2. The second set averages the pressure profiles of 64 mas- 
sive clusters with M > 3 x 1O 14 M from a simulation by 
Dolaget al. (2012). 

3. The third comprises the 40 most massive clusters (i.e., M > 
3 x 1O 14 M ) provided from the numerical simulations de- 
scribed in Battaglia et al. (2010, 201 1). 

All simulations include treatment of radiative cooling, star 
formation and energy feedback from supernova explosions. 
Simulations by Dolag et al. (2012) and Battaglia et al. (2010, 
2011) have prescriptions for AGN feedback. The different sim- 
ulation sets use different techniques and different implementa- 
tions of the physical processes. This ensures a fair description 
of up-to-date theoretical predictions, hopefully bracketing the 
plausible range on the thermal pressure profiles distribution in 
clusters of galaxy. 

We kept the B04+N07+P08 sample as used by A 10, even 
though its mass limit is slightly smaller than our ESZ-XMM 
sample. The higher mass threshold for the two other samples 
enforce the best possible match to our data. As for the observed 
profiles, the simulated cluster profiles were renormalised by the 
factor f(M). The average value for this factor is 1 .0, 1.1 and 1 .05 
for the B04+N07+P08, Dolag et al. (2012) and Battaglia et al. 
(2010) simulations, respectively. Finally, we accounted for the 
differences in the definition of R^m and M500 between obser- 
vations and simulations. We corrected from the differences be- 
tween the true mass and the hydrostatic masses which impact the 
estimation of both Z?5oo and P500 (see Eq. 10). We refer to Sec 4. 1 
of A10 for a complete discussion. Assuming a hydrostatic bias 
r cor = M^/Mire = 1-15 (Kay 2004; Piffaretti & Valdarnini 
2008), we divided the values of R500 and P500 derived from the 
simulations by r\£ and r^ r 3 , respectively. 

All three sets of data are reported in Fig. 6. It is well be- 
yond the scope of this paper to extensively discuss the compar- 
ison between theoretical predictions. Nevertheless, we note that 
they agree within their respective dispersions across the whole 
radial range. The Dolag et al. (2012) and Battaglia et al. (2010, 
2011) profiles best agree within the central part, and are flatter 
than the B04+N07+P08 profile. This is likely due to the imple- 
mentation of AGN feedback, which triggers energy injection at 
the clusters centre, balancing radiative cooling and thus stopping 
the gas cooling. In the outer parts where cooling is negligible, 
the B04+N07+P08 and Dolag et al. (2012) profiles are in perfect 
agreement. The Battaglia et al. (2010) profile is slightly higher, 
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Fig. 6: Measured pressure profiles from Planck and XMM-Newton, dis- 
played as in Fig. 4. The three shaded areas mark the dispersions about 
the average profiles for three samples of simulated clusters: in blue, 
the simulations from Borgani et al. (2004), Nagai et al. (2007a) and 
Piffaretti & Valdarnini (2008), which were used in A10 to derive the 
universal pressure profile together with the REXCESS data; in green, 
the simulated sample of clusters from Dolag et al. (2012); and in brown 
the simulated clusters from Battaglia et al. (201 1). The corresponding 
average profiles for each set of simulations are plotted as solid lines 
using the same colour scheme. In the bottom panel, we present the x 
profiles between the observed profile and the simulated average pro- 
files, taking into account their associated dispersion. See Sect. 6.2 for 
more details. 



but still compatible within its dispersion with the two other sets. 
Here again it might be due to the specific implementation of the 
simulatons. 

7.2.2. Comparison with the Planck pressure profile 

Figure 6 compares our SZ+X-ray observational constraints to 
the theoretical predictions described above. Overall, our ob- 
served pressure profile lies within the scatter across the profiles 
from these various sets of simulated clusters. 

In the central parts, the data points lie on the lower envelope 
of the scatter of the B04+N07+P08 simulations, similarly to the 
REXCESS sample profiles (left panel of Fig. 4 and figure 7 in 
A10). This flatter inner slope is more compatible with that of 
the clusters from simulations which include the effect of AGN 
feedback (Dolag et al. 2012; Battaglia et al. 201 1), although our 
measured points lie above them. 

With increasing radius (i.e., R > R^oo) both SZ and X- 
ray profiles are marginally compatible with the stacked pro- 
files for simulated clusters. Interestingly our profile lies above 
the B04+N07+P08 and Dolag et al. (2012) simulated clusters, 
and below the Battaglia et al. (201 1) clusters. In the outer parts 
of halos, the predicted pressure in numerical simulations is es- 
sentially sustained by gravitational heating. The general agree- 
ment (within the dispersion) with our measured pressure pro- 
file favours the equilibrium between ions and electrons in the 
outskirts, i.e., both populations have an equivalent temperature. 



Indeed, a difference would result in a drop in thermal pressure 
(Rudd & Nagai 2009). This may suggest that we have the cor- 
rect global understanding of the intra-cluster thermal pressure 
support. However, the discrepancies between simulations and 
the average observed profile in the clusters' outer parts calls for 
a deeper investigation of the gas distribution and structure in 
clusters. Thus we need more detailed modelling of the baryon 
physics in cluster outskirts, where incomplete virialisation of 
the intra-cluster medium affects the thermal pressure. Even if a 
universal behaviour of cluster quantities is anticipated from the 
simplest model of gravitational collapse, dispersion is also ex- 
pected, even without considering the effects of non-gravitational 
physics, as shown by numerical simulations of halo formation 
(e.g., Navarro et al. 1997; Reed et al. 201 1; Bhattacharya et al. 
201 1; Gao et al. 2012). This is plainly reflected in the large level 
of dispersion that we observe across our sample of simulated 
clusters. 

On the observational side, our investigation of the sub- 
sample of CC and non-CC clusters did not show differences in 
the outer shape of the pressure profile (CC and non-CC clus- 
ter profiles are compatible within the measured statistical error 
bars). This might indicate that the aforementioned dispersion 
is not mainly driven by the dynamical state of clusters, which 
likely affects primarily the inner parts of halos. If our Planck 
constraints over the whole ESZ-XMM sample give a flatter outer 
shape with respect to the average simulated pressure profiles by 
B04+N07+P08 and Dolag et al. (2012), we need nonetheless to 
keep in mind that our current sample may be affected by selec- 
tion biases, as discussed in Planck Collaboration XI (201 1). 

7.3. Constraint on the gas mass fraction 

As already mentioned, X-ray measurements hardly reach den- 
sity contrasts of ~ 200 (e.g., George et al. 2009; Reiprich et al. 
2008; Urban et al. 201 1; Simionescu et al. 201 1), so modelling 
beyond R200 only relies on predictions from numerical simula- 
tions. We used our Planck and XMM-Newton derived pressure 
profile to investigate the average gas mass fraction distribution 
across our sample. Assuming the ICM to be a perfect gas, the 
thermal pressure is the genuine product of the temperature and 
the density, i.e., P(r) cc n e (r) x kT(r). Following previous X- 
ray works (Pratt et al. 2007; Arnaud et al. 2010; Democles et al. 
2010), we used the best fit to the spectroscopic temperature pro- 
file (Democles et al. 2012) with the analytical function proposed 
by Vikhlinin et al. (2006) for each of our clusters. Beyond the 
reach of the X-ray observations (i.e., [0.5-1.2] R500), we have 
worked with two hypotheses for the temperature: (1) the extrap- 
olation of the best fit model; or (2) a constant temperature fixed 
to the last radial bin value. We have then derived the density pro- 
files, and integrated them over the cluster volume to obtain the 
gas mass radial distribution. Meanwhile, we modelled each total 
mass profile with an NFW profile with the mass, M500, derived 
from Eq. 9, and the concentration estimated from the c — M rela- 
tion by Bhattacharya et al. (201 1). We checked that the average 
of the individual NFW mass models agrees with the NFW model 
derived from the average values of M500 and C500 across the sam- 
ple. The average gas fraction profile is computed as the ratio of 
the average gas mass model and the total mass model profiles. 

Figure 7 shows the resulting gas mass fraction distribution 
for the two hypotheses on the temperature at large radii (red and 
blue solid lines). The red and blue striped areas picture the sta- 
tistical errors measured from the SZ+X-ray pressure profile. The 
red and blue shaded areas overplot the XMM-Newton and Planck 
dispersion onto the pressure profile (see Fig. 4) using hypothesis 
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Fig. 7: Gas mass fraction profile derived from the combined Planck 
and XMM-Newton pressure profile, assuming for the temperature pro- 
file: (1) the average best fit model across the sample from X-ray spec- 
troscopy (red line and striped area); or (2) the same but extrapolat- 
ing beyond ^ 50 o to a constant value equal to the average tempera- 
ture measured in the last radial bin across the sample (blue line and 
striped area). The green dashed curve marks the expected gas frac- 
tion profile assuming the A10 pressure profile and kT(r) as in hypoth- 
esis (1). The star gives f ais (r = /?soo) for REXCESS clusters with 
M 5 oo > 5 x 10 14 M (Pratt et al. 2009). In maroon we reproduce the 
gas mass fraction profile derived from Suzaku measurements for the 
Perseus cluster (Simionescu et al. 2011). The solid and dashed black 
lines mark the cosmic baryon fraction expected from CMB measure- 
ments (Komatsu et al. 2011) and the expected gas fraction, assuming 
that 12% of baryons are in stars, respectively. The shaded blue and red 
areas translate the dispersion in the pressure profile across the ESZ- 
XMM sample as shown on Fig. 4, assuming hypothesis (1) for kT(r). 



( 1 ) for the temperature, and provide an illustration of the cluster- 
to-cluster variation of / gas within our sample. 

From the stacked extrapolated temperature models, the av- 
erage temperature is 1.6 keV within (2-3) R500, which leads 
to /gas ~ 0.2. An overestimation of this temperature by a third 
(i.e., kT ~ 1 keV) would increase the gas fraction likewise (i.e., 
/gas ~ 0.3). For the lower bound, hypothesis (2) fixes the tem- 
perature of the averaged profile beyond R^qq to ~ 4 keV, which 
likely overestimates the true average value of the temperature 
across our sample in the (2-3)/?5oo radial range. Therefore we 
are confident that both hypotheses fairly bracket the range for 
the gas fraction distribution in massive clusters out to 3/?5oo- 
With hypothesis (1) for kT(r), the green dashed line gives the 
gas fraction for the A10 pressure profile. As for our SZ+X-ray 
pressure profile, it leads to / gas values above the expected gas 
fraction out to 3 /?5oo- 

At ^?5oo, the measurement of / gas for massive clusters pre- 
sented by Pratt et al. (2009) (see also Giodini et al. 2009, for 
massive clusters with (4-lO)xlO 14 M ) is lower by a factor 1.3 
with respect to our estimate (although compatible within +2<x). 
We can also compare our result to the gas fraction profile re- 
ported by Simionescu et al. (2011) using Suzaku measurements 
in the direction of the Perseus cluster. Note that this time con- 
suming X-ray observation sampled less than 5% of the sky area 



covered by the cluster. This limited spatial coverage consid- 
ered along with our plausible dispersion envelope (red and blue 
shaded areas on Fig. 7), make this specific measurement neither 
representative nor out of the allowed range for the gas fraction 
in cluster outskirts. 

The constrained interval for the gas fraction from our ob- 
served pressure profile encompasses the cosmic gas fraction de- 
rived from CMB measurements (solid black line on Fig. 7, from 
WMAP-1 results by Komatsu et al. 201 1) and the expected hot 
gas fraction (dashed black line) assuming that 12% of cluster 
baryons are in stars (Gonzalez et al. 2007; Giodini et al. 2009). 



8. Conclusion 

From the Planck nominal mission (i.e., 14 months of survey), 
making use of its full spectral coverage of the SZ spectrum, we 
have extracted and reconstructed the SZ signal distribution in 
clusters of galaxies for a sample of 62 massive nearby clusters. 
All are individually detected in the survey with high significance 
(Planck Collaboration VIII 2011) and were previously used to 
investigate the total integrated SZ flux and the SZ scaling rela- 
tions (Planck Collaboration XI 201 1). We have scaled and aver- 
aged together all the SZ profiles in the sample into a stacked pro- 
file. We have statistically detected the SZ signal out to 3 x 7?goo, 
providing, for the first time, stringent observational constraints 
on the ICM gas beyond R^m and out to a density contrast of 
8~ 50-100. 

From the Planck SZ raw profiles, we have derived the un- 
derlying thermal pressure profiles of the ICM. Together with 
the pressure profile derived from the XMM-Newton data, we 
have provided for the first time a comprehensive observational 
view of the distribution of thermal pressure support in clusters 
from 0.01 out 3 x /?5oo- We have fitted these unique measure- 
ments to a generalised NFW profile. Our best analytical repre- 
sentation over this wide radial range is given by the parameters 
[Po,c 500 ,y,a,/3] = [6.41,1.81,0.31,1.33,4.13]. 

Our observational measurements further confirm the agree- 
ment of the SZ and X-ray constraints on the intra-cluster 
gas properties within the inner part of the clusters com- 
monly considered to be virialised, i.e., R500, as found in 
Planck Collaboration X (2011) and Planck Collaboration XI 
(2011). Overall it also agrees with a wide range of simulated 
clusters through the whole radial range, although in the central 
regions it matches best the numerical predictions that implement 
AGN feedback. The statistical nature of our stacked detection at 
large radii provides the average trend for the thermal pressure 
distribution in cluster outskirts out to 3 x R500, which is slightly 
flatter than most theoretical predictions. 

In conjunction with X-ray constraints on the temperature 
profile of our clusters, we have derived the profile of the gas mass 
fraction out to the cluster outskirts. From reasonable hypothe- 
ses on the gas temperature to account for the lack of constraints 
beyond R500, we have bracketed a range for the gas fraction in 
the cluster outer regions, which is compatible with the cosmic 
baryon fraction and the expected gas fraction in halos. 

The processes governing the thermodynamical state of 
the outer regions in clusters still need to be understood 
from the theoretical and observational points of view (see 
Kravtsov & Borgani 2012, for a review and references). Issues 
such as gas clumping will affect the pressure estimation (i.e., ef- 
fect on the X-ray surface brightness, see Roncarelli et al. 2006), 
departures from hydrostatic equilibrium and contribution from 
non-thermal pressure (due to magnetic fields and/or cosmic rays) 
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will modify the gas fraction. In this regard SZ observations pro- 
vide a straightforward description of the thermal pressure dis- 
tribution in massive halos, and a clear path to the gas fraction 
determination. Therefore with Planck we bring unique observa- 
tional constraints that are extremely valuable to further test and 
understand the physics at play in the outskirts of clusters. 
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Appendix A: SZ map reconstruction methods 

In this Appendix we discuss the different SZ map reconstruction 
methods applied to the Planck data in Sect. 4. 1 . 

1. MILCA (Hurieretal. 2010): The thermal SZ signal re- 
construction is performed on the six Planck all-sky maps 
from 100 GHz to 857 GHz. MILCA (Modified Internal 
Linear Combination Algorithm) is a component separation 
approach aiming at extracting a chosen component (here 
the thermal SZ signal) from a multi-channel set of in- 
put maps. It is mainly based on the well known Internal 
Linear Combination approach (e.g., Eriksen et al. 2004), that 
searches for the linear combination of the input maps that 
minimises the variance of the final reconstructed map by im- 
posing spectral constraints. In this paper, we applied MILCA 
using two constraints: preservation of the thermal SZ sig- 
nal (knowing the SZ spectral signature); and removal of the 
CMB contamination in the final SZ map (also making use of 
the well known spectrum of the CMB). In addition, to com- 
pute the weights of the linear combination, we have used the 
extra degrees of freedom to minimise residuals from other 
components (2 degrees) and from the noise (2 degrees). The 
noise covariance matrix was estimated from the frequency 
error maps (see Sect. 2). 

2. NILC: Needlet ILC performs a linear combination of the 
observed maps which has minimum variance under the con- 
straint of offering unit response to the component of interest 
(here the thermal SZ, whose frequency scaling is known). 
The weights of the ILC depend on the covariance between 
the various observations, and can be computed, for exam- 
ple, on domains of the observed pixels (pixel space) or an- 
gular scales (spherical harmonics). In the case of NILC, co- 
variances (and hence weights for component separation) are 
computed independently in domains of a needlet decom- 
position (spherical wavelet frame). The needlet decompo- 
sition provides localisation of the ILC filters both in pixel 
and in multipole space, allowing us to deal with local con- 
tamination conditions varying both in position and in scale. 
NILC was developed to extract a CMB map from WMAP 
data (Delabrouille et al. 2009) and was also tested for SZ ef- 
fect extraction in Leach et al. (2008). Multi-component ex- 
tensions have been investigated by Remazeilles et al. (201 1). 

3. GMCA (Bobin et al. 2008) is a blind source separation 
method developed for separating sources from instantaneous 
linear mixtures. The components are assumed to be sparsely 
represented (i.e., have a few significant samples in a spe- 
cific basis) in a so-called sparse representation *P (typically 
wavelets). The assumption that the components have a sparse 
representation in the wavelet domain is equivalent to as- 
suming that most components have a certain spatial regular- 
ity. These components and their spectral signatures are then 
recovered by minimising the number of significant coeffi- 
cients in W, Recently, L-GMCA has been further introduced 



to analyse the CMB data in a local and multi-scale manner 
(Bobin et al. 2011). More precisely, the multi-channel data 
are analysed in four frequency bands in spherical harmonics 
(i.e., wavelet bands). In each wavelet band, GMCA is ap- 
plied locally on small patches with band-dependent sizes. In 
each band, the observations are analysed at the same band- 
dependent resolution. In Bobin et al. (2011), it was shown 
this local analysis approach enhances the separation qual- 
ity. The spectral signatures of CMB, free-free and SZ are 
assumed to be known. 



10° r 




R/R500 



Fig.A.l: Comparison between the SZ profile reconstruction from 
MILCA, NILC and GMCA The comparison is presented in the form 
of the averaged stacked radial profile from our 62 clusters of the ESZ- 
XMM sample, individually scaled respectively in the x- and y-axis di- 
rections according to R^oo and C> (Eq 12). The light-red shaded area 
marks the dispersion about the average stacked profile for the MILCA 
method. The points within each individual profile are correlated at about 
the ~ 20% level (see Sect. 4.2) and the plotted errors correspond to the 
square root of the diagonal elements of the covariance matrix of the pro- 
file. The solid black line (labelled in the legend as 'Expected from P A 10 ') 
is the average stacked profile obtained from the expected SZ individ- 
ual profiles drawn from the universal pressure profile by A 10, param- 
eterised according R^oo and F500 derived from the XMM-Newton data 
analysis presented in Planck Collaboration XI (201 1)). The two dotted 
lines indicate the associated dispersion across the sample. The bottom 
panel shows the value of x at each point of the observed profiles with 
respect to the expectations from the universal profile of A10. 

Appendix B: Validation on simulations 

B.1. Profiles 

We have validated our reconstruction of the SZ signal and profile 
computation methods with simulations. We used the Planck sky 
model (PSM, Delabrouille et al. 2012) to simulate the whole 
sky as seen by Planck. The PSM includes most astrophysical 
components acting as foreground or background contamination 
with respect to the SZ signal (i.e., Galactic dust, free-free and 
synchrotron emission, extragalactic sources in the infrared and 
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radio parts of the spectrum and CMB). For the instrumental 
noise we added to the simulation the all-sky error map drawn 
from the Planck jack-knife maps at each frequency (see Sect. 2), 
thus mimicking the noise properties of the actual survey (noise 
and systematic effects). For the SZ component, we built a full 
sky template for all the 1743 clusters from the MCXC meta- 
catalogue (Piffaretti et al. 201 1), which includes all 62 clusters in 
our sample. Each cluster was modelled assuming spherical sym- 
metry and a thermal pressure distribution following the universal 
profile from A 10. We recall that the universal pressure profile 
provided by these authors is the best fit of a generalised NFW 
profile (Nagai et al. 2007a) to the median profile derived from 
the REXCESS pressure profile and the predictions from three 
different sets of numerical simulations (i.e., from Borgani et al. 
2004; Nagai et al. 2007a,b; Piffaretti & Valdarnini 2008). In our 
case, the universal pressure profile is parametrised using the val- 
ues of ^500 and M 500 from the MCXC, together with the M-Y x 
scaling relation from Arnaud et al. (2007, 2010) to estimate i^oo- 
We processed the simulated Planck sky maps through each 
of our SZ map reconstruction methods. For our sample of 62 
clusters, we extracted the patches and computed the individual 
profiles as described in Sect. 4.2. We stacked the simulated pro- 
files for the three methods after renormalisation, making use of 
the aforementioned Rsqq and Fsoo values from the MCXC pre- 
scription. 

The resulting stacked profiles are displayed in Fig. A. 1. All 
three tested methods agree remarkably well over the whole radial 
range. Taking into account the correlated errors of each profile, 
thex 2 over the radial range (0-10)x^ 50 o for NILC and GMCA 
with respect to MILCA are 0.78 and 1.54, respectively. The re- 
constructed profiles are also in good agreement with the recon- 
structed stacked model. Accounting for the errors on the model 
and the correlated errors between points on the reconstructed 
profiles, the reduced^ 2 for MILCA, NILC and GMCA with re- 
spect to the stacked model profile are 0.80, 0.66 and 1.03 within 
3 x R50Q. This excellent agreement also translates into a small 
relative error below 15% in the same range of radii. Beyond 
R ~ 2 R500, the three reconstructed profiles drop slightly below 
the input model, even though they still agree within their dis- 
persions and errors. This drop is caused by an intrinsic over- 
prescription of point sources in the PSM simulations with re- 
spect to real data. Despite our careful masking, the unidentified 
or unresolved extra sources in the PSM affect the power spec- 
trum of the noise, increasing the background level and making it 
more difficult to detect the SZ signal further out in the simula- 
tions than in real data. We stress that this is an effect limited to 
our test simulations, and not relevant for the real data. 

B.2. Consistency between F500 measurements 

To further validate the above, and as a consistency check with 
previous Planck results, we fitted each individual SZ profile with 
a projected, PSF-convolved universal pressure profile. We fixed 
/?50o to the best fitting X-ray value from Planck Collaboration XI 
(201 1) and only fitted the normalisation, Ysoq. For a given pro- 
file, P, with an associated covariance matrix, C, and a model M, 
thex 2 statistic can be expressed as 

X 1 = (P ~ F500 x M) T C-\P - F 500 x M). (B. 1) 

The solution that minimises x 2 , isoo> and its associated uncer- 
tainty, <Ty 500 , are analytic: 

F 5 oo = oi oo M T C- 1 P; (B.2) 
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Fig.B.l: Comparison of Y 5m values from the Planck cluster detec- 
tion method (MMF3) against our best fit for each of our 62 clusters 
(from NILC in this case). In both cases the distribution of thermal pres- 
sure is assumed to follow a universal pressure profile (A10) for which 
the scaled radius R 5m is fixed to the XMM-Newton best fit value. See 
Sect. 6.3 for details. 



<r 2 5oo =(M T C" 1 M)- 1 . (B.3) 

We compared the resulting value of F500 to those obtained 
from the detection algorithms PwS and MMF3. Both detec- 
tion algorithms were run under the same assumptions, i.e., uni- 
versal pressure profile with the "non-standard" parameterisa- 
tion and a fixed size (Rsoo). The fluxes were derived from the 
Planck nominal mission survey maps. These fluxes are in agree- 
ment with those extracted from the Planck first year survey 
(Planck Collaboration XI 2011). For instance, for MMF3, the 
average ratio is 1 .04 + 0.09, compatible within the average rela- 
tive errors with the fluxes extracted from the nominal mission. 

The relation between the integrated Comptonisation param- 
eters derived from the SZ profiles and from the detection meth- 
ods is displayed for MMF3 vs NILC fuxes in Fig. B.l. A BCES 
linear fit accounting for the errors in both x- and y-axis direc- 
tions (Akritas & Bershady 1996) yields a best fitting slope and 
normalisation of 1.00 + 0.02 and 0.0 1+0.01 0j> 2 = 0.83), respec- 
tively. There is 4% dispersion about this relation. Furthermore, 
the median and associated deviation for the one-to-one ratios of 
PwS and MMF3 fluxes compared to the fitted NILC profile val- 
ues are 1.03 + 0.11 and 0.98 + 0.08, respectively. The median 
ratios for the errors are 1.10 + 0.36 are 1.64 + 0.34. The errors as- 
sociated with our best fit values are in very good agreement with 
those from PwS; however, they are smaller than those derived 
from the MMF3 method. This difference with MMF3 might be 
explained by the different treatment of the noise. Whilst we es- 
timated the noise from the surroundings of the clusters, defined 
according to the scaled radius of each target, and PwS also used 
a local computation of the noise, the MMF3 algorithm used a 
larger area to characterise the noise properties. Our tests against 
the PwS and MMF3 fluxes emphasise our grasp on the charac- 
terisation of the local noise for each target, and further assess the 
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reliability of our estimates of the statistical errors in our stacked 
profiles. 

Such agreement between the fluxes derived from our fit 
over the cluster profiles and those from the detection algo- 
rithm, strongly confirms the reliability of our SZ flux estimate 
and further strengthen the overall result stated in the Planck 
collaboration early papers that the SZ and the X-ray fluxes 
within /?5oo are fully consistent (Planck Collaboration VIII 201 1; 
Planck Collaboration IX 2011; Planck Collaboration X 2011; 
Planck Collaboration XI 201 1). 

Appendix C: Pressure profiles of individual clusters 

In this appendix we present the individual pressure profiles ob- 
tained from Planck and XMM-Newton data for each of the 62 
clusters in our sample (see Fig. C.l). We show the agreement 
with the universal pressure profile (A10) and provide the best 
fitting model for each cluster to the GNFW profile in Table C. 1 . 
The following hard limit priors were adopted during the fitting 
procedure: < P () < 100 ; < C500 < 10 ; < a < 10 ; and 
<J3 < 15. 
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Table C.l: Best fit parameters for a generalised NFW pressure profile for individual clusters. 



Cluster name 


Pq 


£'500 


r 


a 




A2163 


5.28 


3.64 


0.31 


2.51 


2.87 


RXCJ0014.3-3022 


2.74 


1.05 


0.31 


1.83 


7.24 


A2204 


35.29 


2.79 


0.31 


0.83 


3.68 


RXJ1720.1+2638 


30.51 


2.15 


0.31 


0.74 


3.72 


A2034 


4.92 


3.82 


0.31 


10.00 


2.62 


A2261 


22.88 


4.97 


0.31 


1.25 


2.79 


MACSJ2243. 3-0935 


3.41 


2.53 


0.31 


4.27 


3.38 


A1914 


15.57 


5.83 


0.31 


2.75 


2.66 


A2390 


5.96 


0.26 


0.31 


0.90 


14.94 


A1763 


5.61 


1.43 


0.31 


1.10 


4.00 


A2218 


9.65 


5.51 


0.31 


2.99 


2.23 


RXCJ0043.4-2037 


6.72 


2.92 


0.31 


2.15 


3.52 


A1576 


9.23 


2.81 


0.31 


1.45 


3.45 


A665 


3.13 


0.14 


0.31 


0.80 


14.38 


A773 


7.62 


2.60 


0.31 


1.46 


3.29 


A963 


8.82 


0.17 


0.31 


0.71 


13.04 


A781 


1.78 


1.82 


0.31 


5.56 


3.77 


A520 


3.31 


2.75 


0.31 


4.00 


2.98 


A1413 


17.94 


1.62 


0.31 


0.83 


4.31 


RXCJ0532.9-3701 


11.83 


3.81 


0.31 


2.00 


3.39 


RXCJ0528.9-3927 


4.62 


0.07 


0.31 


0.64 


14.74 


A868 


5.65 


1.88 


0.31 


1.48 


3.94 


RXCJ0232.2-4420 


11.50 


1.95 


0.31 


1.10 


4.26 


RXCJ05 16.7-5430 


1.36 


0.24 


0.31 


1.21 


14.75 


RXJ0658-55 


5.69 


2.92 


0.31 


2.91 


3.46 


RXCJ1131.9-1955 


4.18 


2.03 


0.31 


1.90 


3.85 


RXCJ0303. 8-7752 


4.23 


3.24 


0.31 


4.94 


2.97 


A1650 


14.21 


0.78 


0.31 


0.78 


6.17 


A1651 


15.87 


3.89 


0.31 


1.23 


2.81 


A1689 


33.95 


1.76 


0.31 


0.77 


4.49 


A3921 


6.09 


0.75 


0.31 


0.96 


5.63 


A3827 


9.74 


3.42 


0.31 


1.54 


2.97 


A3911 


4.39 


2.30 


0.31 


1.98 


3.07 


AS 1063 


11.80 


1.30 


0.31 


1.08 


6.18 


A2219 


7.04 


3.25 


0.31 


1.89 


2.90 


A2255 


1.82 


0.53 


0.31 


1.41 


8.35 


R0605.3518 


11.25 


0.07 


0.31 


0.58 


14.91 


R2218.3853 


7.51 


1.97 


0.31 


1.35 


4.21 


R0547.3152 


7.13 


2.22 


0.31 


1.77 


4.59 


R2217.3543 


13.20 


2.30 


0.31 


1.12 


3.97 


R2234.3744b 


9.98 


4.10 


0.31 


2.66 


3.16 


R0645.5413c 


6.54 


0.49 


0.31 


0.88 


8.03 


RXJ2228 


6.92 


0.92 


0.31 


0.99 


5.85 


RXJ1206 


6.29 


0.12 


0.31 


0.70 


14.80 


A401 


5.80 


3.79 


0.31 


2.08 


2.46 


A478 


30.40 


3.00 


0.31 


0.84 


3.53 


A2065 


3.73 


0.35 


0.31 


1.09 


15.00 


A1795 


8.11 


0.10 


0.31 


0.63 


14.98 


A2029 


21.48 


0.91 


0.31 


0.66 


5.29 


ZwC11215 


4.88 


2.46 


0.31 


1.65 


3.17 


A85 


5.99 


0.02 


0.31 


0.48 


14.97 


A3558 


6.04 


1.77 


0.31 


1.12 


3.58 


A2256 


2.72 


1.65 


0.31 


2.41 


4.38 


A3112 


24.16 


0.03 


0.31 


0.44 


14.11 


A3158 


5.93 


1.63 


0.31 


1.17 


4.11 


A1644 


2.08 


0.03 


0.31 


0.60 


14.89 


A754 


1.76 


2.42 


0.31 


2.63 


3.66 


A3266 


3.05 


1.15 


0.31 


1.55 


5.60 


A3376 


1.49 


1.42 


0.31 


3.57 


4.89 


A3528s 


5.72 


0.01 


0.31 


0.36 


11.63 


A3532 


4.79 


1.90 


0.31 


1.08 


2.94 


A0119 


2.38 


1.67 


0.31 


1.81 


3.44 
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(a) 

Fig. C.l: Measured individual pressure profiles for the ESZ-XMM sample using the Planck and XMM-Newton data. As for the stacked pressure 
profile, the maximum difference point-to-point between the three SZ reconstruction methods have been taken into account in the Planck measure- 
ment error budget (Sect. 6.3). The best fit model on each individual profile is shown as a blue solid line. The black and green solid curves mark 
the best fit model to the stacked SZ+X-ray pressure profile (Sect. 6.2) and the A10 best fit model, respectively 
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Fig.C.l: continued. 
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Fig.C.l: continued. 
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Fig.C.l: continued. 



